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I. INTRODUCTION 



In the past few years there has been considerable work on non-Hermitian Hamiltonians. 
Among this kind of Hamiltonians, much attention has been focused on the investigation of 
properties of so-called PT-symmetric Hamiltonians. Following the early studies of Bender 
et al. [1], the 'PT-symmetry formulation has been successfuly utilized by many authors 
[2-8]. The PT-symmetric but non-Hermitian Hamiltonians have real spectra whether the 
Hamiltonians are Hermitian or not. Non-Hermitian Hamiltonians with real or complex 
spectra have also been analyzed by using different methods [3-6,9]. Non-Hermitian but VT- 
symmetric models have applications in different fields, such as optics [10], nuclear physics 
[11], condensed matter [12], quantum field theory [13] and population biology [14]. 

Exact solution of Schrodinger equation for central potentials has generated much interest 
in recent years. So far, some of these potentials are the parabolic type potential [15], the 
Eckart potential [16,17], the Fermi-step potential [16,17], the Rosen- Morse potential [18], 
the Ginocchio barrier [19], the Scarf barriers [20], the Morse potential [21] and a potential 
which interpolates between Morse and Eckart barriers [22]. Many authors have studied on 
exponential type potentials [23-26] and quasi exactly solvable quadratic potentials [27-29]. 
In addition, Schrodinger, Dirac, Klein-Gordon, and Dufifin-Kemmer-Petiau equations for a 
Coulomb type potential are solved by using different method [30-34]. The exact solutions 
for these models have been obtained analytically. 

Further, using the quantization of the boundary condition of the states at the origin, 
Znojil [35] studied another generalized Hulthen and other exponential potentials in non- 
relativistic and relativistic regions. Domingues-Adame [36] and Chetouani et al. [37] also 
studied relativistic bound states of the standard Hulthen potential. On the other hand, Rao 
and Kagali [38] investigated the relativistic bound states of the exponential-type screened 
Coulomb potential by means of the one-dimensional (1-D) Klein-Gordon equation. How- 
ever, it is well-known that for the exponential-type screened Coulomb potential there is no 
explicit form of the energy expression of bound-states for Schrodinger [39], KG [38] and 
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also Dirac [16] equations. §im§ek and Egrifes [31] have presented the bound-state solutions 
of ID Klein-Gordon (KG) equation for PT-symmetric potentials with complex generalized 
Hulthen potential. In a latter study, Egrifes and Sever [32] investigated the bound-state 
solutions of the ID Dirac equation with PT-symmetric and non-PT-symmetric real and 
complex forms of the generalized Hulthen potential. Yi et al. [40] obtained the energy 
equations in the KG theory with equally mixed vector and scalar Rosen-Morse-type poten- 
tials. Berkdemir et al. [41] obtained the bound-state eigenvalues and eigenfunctions for the 
Schrodinger equation using a generalized form of Woods-Saxon (WS) potential by means of 
Nikoforov-Uvarov method [42] . In recent works, we have investigated the £ ^ bound-state 
solutions of the ID Schrodinger equation with the real and complex forms of the modified 
Hulthen and WS potentials [43,44] for their energy spectra and correponding wave functions 
using the Nikiforov-Uvarov (NU) method. In addition, we have extended our study to rela- 
tivistic models in solving the spinless Salpeter equation analytically for its exact bound-state 
spectra and wavefunctions for the real and complex forms of the "PT-symmetric generalized 
Hulthen potential [45]. 

In the present work, we investigate the bound-state solutions of the ID Klein-Gordon 
(KG) equation with real and complex forms of the generalized Woods-Saxon (WS) potential. 
Hence, the objective of our work is to determine the exact energy levels of relativistic KG 
particles trapped in a spherically symmetric generalized WS potential which possesses a 
PT-symmetry as well. In this context, it is possible to convert KG to a Schrodinger-like 
equation which takes the hypergeometric form to which one may apply the NU method. 

The organization of this work is as follows. After a brief introductory discussion of the 
NU method in Section II, we discuss the KG problem and then obtain the exact bound-state 
energy spectra for real and complex cases of generalized WS potentials and their correspond- 
ing eigenfunctions in Section III. In Section IV, we discuss our results and mention briefly 
the scope of this relativistic study. 
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II. THE NIKIFOROV-UVAROV METHOD 



In this section we outline the basic formulations of the method. The Schrodinger equation 
and other Schrodinger-type equations can be solved by using the Nikiforov-Uvarov (NU) 
method which is based on the solutions of general second-order linear differential equation 
with special orthogonal functions [42]. It is well known that for any given ID radial potential, 
the Schrodinger equation can be reduced to a generalized equation of hypergeometric type 
with an appropriate transformation and it can be written in the following form 

a(s) a z {s) 

where a(s) and a(s) are polynomials, at most of second- degree, and f(s) is of a first-degree 
polynomial. To find particular solution of Eq.(l) we apply the method of separation of 
variables using the transformation 

^n(s) = (f> n (s)y n (s), (2) 

which reduces Eq.(l) into a hypergeometric-type equation 

cr(s)y:(s) + r(s)y' n (s) + \y n (s)=0, (3) 

whose polynomial solutions y n (s) of the hypergeometric type function are given by Rodrigues 
relation 

Vn{s) = ^y^T > (n = Q, 1, 2, ...) (4) 

where B n is a normalizing constant and p{s) is the weight function satisfying the condition 
[42] 

where w(s) = a(s)p(s). On the other hand, the function 0(s) satisfies the condition 



where the linear polynomial 7r(s) is given by 



^^Mj 2 _ 5(s) + Ms)> (7) 



from which the root is the essential point in the calculation of 7r(s) is determined. Further, 
the parameter A required for this method is defined as 

\ = k + 7i'{s). (8) 

Further, in order to find the value of k, the discriminant under the square root is being set 
equal to zero and the resulting second-order polynomial has to be solved for its roots 
Thus, a new eigenvalue equation for the SE becomes 

X n + nr'(s) + n{n ~ l) a"{s) =0, (n = 0, 1, 2, ...) (9) 

where 

r(s) =f(s) + 27r(s), (10) 
and it must have a negative derivative. 

III. EXACT BOUND STATE SOLUTIONS OF THE GENERALIZED 
WOODS-SAXON POTENTIAL 

In this section, we formulate the Nikiforov-Uvarov solution for the relativistic motion of 
a spin-zero particle bound in spherically symmetric well-known spherical symmetric gener- 
alized Woods-Saxon potential of the form [46] 

V q (r, Ro) = -V 6 ° v q > 0, R > a, (11) 
1 + qe i^^> 

where r refers to the center-of-mass distance between the projectile and the target nuclei 
(the r range is from to oo). The relevant parameters of the nuclear potential are given as 
follows: Rq = t-qA 1 / 3 is to define the confinement barrier position value of the corresponding 



spherical nucleus or the width of the potential, A is the target mass number, r is the 
radius parameter, V controls the barrier height of the Coulombic part, a is the surface 
diffuseness parameter has to control its slope, which is usually adjusted to the experimental 
values of ionization energies. Further, q is a shape (deformation) parameter, the strength 
of the exponential part other than unity, set to determine the shape of potential and is 
arbitrarily taken to be a real constant within the potential. Further, we remark that the 
spatial coordinates in the potential are not deformed and thus the potential still remains 
spherical. Obviously, for some specific q values this potential reduces to the well-known 
types, such as for q = to the exponential potential and for q = —1 and a = S^ 1 to the 
modified Hulthen potential (cf. [43,44,45] and the references therein). 

For a scalar particle of rest mass m and total energy E nh the radial part of the KG 
equationin in three-dimensional spherical coordinates [47] 



h 2 d 2 h 2 i(i + i) 



+ 



+ 



2 4 

m c 



[Era - V{r)\ 



i>ni{r) = 0, 



(12) 



2m dr 2 2mr 2 2mc 2 
where if) n i(r) is the reduced radial wave function. 

To calculate the energy eigenvalues and the corresponding eigenf unctions, we substitute 
the Hermitian real- valued ID generalized WS potential (the r range is from to oo): 



V q (x) 



■Vn 



1 + qe~ 



-, a = 1/a, x = r — R 



(13) 



into the ID VT— symmetrical Hermitian spinless KG equation (12) for / = case (i.e., 
s-wave states) and then obtain 



d 2 A iq (x) 1 



dx 2 



h 2 c 2 



[E 2 -m 2 c*) + V 2 - 



-2ax 



+ 2V E n 



i> nq {x) = 0. (14) 



l + qe- ax ) 2 {l + qe- a 

We employ the following novel dimensionless transformation parameter, s(x) = (l + ge~ aa: ) _1 
, which is in real phase and found to maintain the transformed wavefunctions finite on 
the boundary conditions (i.e., — oo <x<oo^0<s<l) [44]. Hence, applying such 
transformation to the upper equation after setting H — c — 1, it gives 



d 2 ip nq (s) 1 - 2sdip ng (s) 



~d? 



+ 



S — S z 



nq \ 

ds 



+ 



a 2 (s — s 2 ) 2 



{E 2 n - m 2 ) + V 2 {1 - s) 2 + 2E n V Q {\ - s) ^ nq {s) = 0, 



(15) 



where V = V /q. With the dimensionless definitions given by 



— e 



[El - m 2 ) 



cr 



>o, p 2 



2E n V 



, i 2 %.(!■:;< »r- >' 2 > o. 7 2 >o), (16) 



cr 



cr 



(17) 



for bound-states, one can arrive at the simple hypergeometric equation given by 

(s — s J (s — s^) 

Hence, comparing the last equation with the generalized hypergeometric type, Eq.(l), we 
obtain the associated polynomials as 

f(s) = 1 - 2s, a(s) = (s - s 2 ), a(s) = s 2 7 2 - s (f3 2 + 2 7 2 ) + (3 2 + 7 2 - e 2 . (18) 



When these polynomials are substituted into Eq.(7), with cr'(s) = 1 — 2s, we obtain 

tt(s) = ± v / -s 2 ( 7 2 + k) + s((3 2 + 2 7 2 + k) + e 2 - /3 2 - 7 2 . 



(19) 



Further, the discriminant of the upper expression under the square root has to be set equal 
to zero. Therefore, it becomes 



fi 2 + 2 1 2 + k -AiY + kW + Y-e 2 ) = 0. 



(20) 



Solving Eq.(20) for the constant k, we obtain the double roots as k\^ = f3 2 — 2e 2 ±2e6, where 
b = \fe 2 — (3 2 — 7 2 = | \/l ~~ 4 7 2 — (2n + 1) — e with n — 0, 1, • • • . Thus, substituting these 
values for each into Eq.(19), we obtain 



7r(s) = ± < 



(b - e) s - b; for fci = /3 2 - 2e 2 + 2eb, 
(b + e)s-b; for k 2 = 1 - 2e 2 - 2eb. 

Hence, making the following choice for the polynomial ir(s) as 

tt(s) = -(6 + e)s + 6, 
for k 2 = f3 2 — 2e 2 — 2e6, giving the function: 



(21) 



(22) 



r(s) = -2(1 + 6 + e)s + 1 + 26, 



(23) 



which has a negative derivative of the form t(s) = —2(1 + b + e) =2n— 1 — y/\ — A^y 2 . 
Thus, from Eqs.(8)-(9) and Eqs.(22)-(23), we find 

A = - 7 2 -(6 + e)(6 + e + l), (24) 

and 

\ n = n 2 + n + 2n (e + b) . (25) 

Therefore, after setting \ n = A and solving for E nq , we find the KG exact binding energy 
spectra as 

Enq = -^ q ± ^WTe'W' (26) 

where 

£ = VV« 2 - W 2 - qa(2n + 1), (27) 

providing that the condition q 2 a 2 > 4Vq must be fulfilled for any possible bound-states. In 
this context, it is worthwhile to point out that §im§ek and Egrifes [31] have recently obtained 
a similar expression for the Hulthen potential. We find the corresponding wavefunctions by 
applying the NU method to find the hypergeometric function y n (s) which is the polynomial 
solution of hypergeometric-type equation (3) described with the weight function [42]. By 
substituting ir(s) and a(s) in Eq.(6) and then solving the first-order differential equation, 
we find 

Ms) = s\l-sY. (28) 

To find the function y nq (s), which is the polynomial solution of hypergeometric-type equa- 
tion, we multiply Eq.(3) by p(s) so that it can be written in self- adjoint form [42] 

(a(s)p(s)y' nq (s))' + Xp(s)y nq (s) = 0, (29) 

where p(s) satisfies the differential equation (cr(s)p(s))' = r(s)p(s) which gives 
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p(s) = s 2b (l - s) 2e . 



(30) 



We then obtain the eigenfunction of hypergeometric-type equation from the Rodrigues re- 
lation given by Eq.(4) as 



y nq (s)=D nq s- 2b (l-s)-*— 



,n+2f> , 



1-S 



,n+2c 



(31) 



where D nq is a normalizing constant. In the limit q — > 1, the polynomial solutions of y n (s) 
are expressed in terms of Jacobi Polynomials, which is one of the classical orthogonal polyno- 
mials, with weight function (30) in the closed interval [0, 1] , giving y n ,\{s) ^ P^ b,2e \l — 2s) 
[42]. The radial wave function ^ nq (s) is obtained from the Jacobi polynomials in Eq.(31) 
and <f>(s) in Eq.(28) for the s-wave functions could be determined as 

iM*) = N nq s-\l - s)-^ [s n+2b (1 - s) M ] = N nq s b (l - S yP^\\ - 2s), (32) 



where b = ^ — ^ m ~ En ; s — (\ -f qe~ ax )~ l and N nq is a new normalization constant. 

We now make use of the fact that the Jacobi polynomials can be explicitly written in 
two different ways [48]: 



(33) 



P M {z) 



T(n + p+l) A (n\T(n + p + v + r + 1) (z- 1 
v + 1) 



(34) 



nW(n + p + v + l)^ Q \rJ T(r + p+l) V 2 

where (™) = r!( "_l r) , = r(r+i)r^r+i) • Usin S Eqs.(33)-(34), we obtain the explicit expressions 
for P^\l -2s) 

Pi 2b ' 2e \l - 2s) = (-l) n r(n + 2b + l)r(n + 2e + 1) 



' —\yq n ~p 



x V" 

P ! (« - P) !r (P + 2e + !) r ( n + 26 - p + 1) 



s n " p (l - s) p , 



p(2 & ,2e) ri _ 2 x ^ + 26+1) " (-l)r g r r(w + 26 + 2e + r+1) 

1 ' r(n + 26 + 2e + l) r\(n - r)\T(2b + r + 1) 



(35) 



(36) 
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, = ^2 ( r r(n + 26 + l)r(n + 2&+l) 2 f " (-l)V^ 



r(n + 2e + 2b + 1) p!(n - p)\T(p + 2e + l)r(n + 2b - p + 1) 

f » (-i)yr(n + 26 + 26 + r + l) | 

l^o r!(n-r)!r(26 + r + l) ^ 

where 

i 

A*(p, r) = y s "+ 26 +^ (1 - s )P+ 2e <fe. (38) 
o 

Using the following integral representation of the hypergeometric function [49] 

fs°°-\l - s)— (1 - qs^ds = 2 F l(ao ,/3o : 7o; g ) g^^^ , 

[Re( 7o ) > Re(a ) > 0, |arg(l - q)\ < tt] (39) 

which gives 

i 

2 F 1 (a ,A):ao + l;g)/ao = / s a °-\l - qs^ds, (40) 

o 

p / a \ r(7p)r(7o - q - /3q) 

r(7o - ao)r(7o -AO 

[Re( 7o - a - AO > 0, Re( 7o ) > Re(A0 > 0, (41) 

for g = 1. Setting ceo = n + 2b + r — p + 1, A) — —p — 2e, and 70 = «o + 1, one gets 

, , x _ 2^0, A) : 7o;g) _ (n + 2b + r-p + l)\(p + 2e)\ 

nq{P,T> a (n + 2b + r-p+l)(n + 2e + r + 2b+l)\ 1 } 

A. Real Potentials 

Consider the parameters Vb, 9, and a given in Eq.(13) are all real, then 
(i) For any given a the spectrum consists of real eigenstate spectra E n (Vo, q, a) depending 
on q. The sign of Vq does not affect the bound states. For positive shape parameter q, it 
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is clear that while V — > 0, E n = my 1 — (f^J tends to the value m for the ground state 
(i.e., n — 0) and 0.866m for the first excited state (i.e., n — 1), etc. In these calculations, 
we have used A c = — = — = — which is the compton wavelength of the KG particles. 

L mc ma jr o jr 

(ii) There exist bound states (real solution) in case if the condition 4Vq < q 2 a 2 is 
achieved, otherwise there are no bound-states. 

(iii) There exist bound states in case if the condition AVq + £ 2 < 16g 2 m 2 is achieved, 
otherwise there are no bound-states. 

Moreover, this condition which gives the critical coupling value turns to be 



i.e., there are only finitely many eigenstates. In order that at least one level might exist 
(n — 0), its necessary that the inequality 



Under V, the spatial coordinates (x, y, z) are replaced by (—x, —y, —z) but r is replaced 
by r and not — r, in the radial wave equation (4). Thus, the s-wave differential equation is 
not VT— symmetric. The radial Schrodinger wave equation becomes a different differential 
equation under the action of the VT— operator and does not go into itself. This means that 
we must solve the problem in ID, on the full plane, say x-direction and not in the radial 
direction r. 




(43) 




(44) 



is fulfilled 



B. Complex Potentials 



1. Non-Hermitian VT Symmetric Generalized Woods-Saxon Potential 



Let us consider the case where at least one of the potential parameters be complex: 
If a is a complex parameter (a — > ia), the potential (13) becomes 
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V q (x) = — — — — — — — — [q + cos(arr) - i ffln(oa;)] = V*(-x), (45) 

q 2 + 2qcos(ax) + 1 q 

which is a PT-symmetric but non-Hermitian. It has real spectrum given by 



E m = ~ + ( vVa 2 + W$ - aq(2n + V 



i 



1 m 



2 



16q2 4V 2 - ^q 2 a 2 + 4V 2 - qa{2n + l))'' 

(46) 



2 



if and only if 16q 2 m 2 < - ^jq 2 a 2 + 4\/ Q 2 - qa{2n + 1)J which gives ^g^y + 2(2^ + 
f (2n+ 1) < + 4V 2 . 

Figs. 1(a) and (b) show the variation of the ground-state (i.e., n = 0) as a function 
of the coupling constant Vq for different positive and negative q, and a = A c . Obviously, 
in Fig. 1(a), the non-Hermitian PT-symmetric generalized WS potential generates real 
and negative bound-states for q > 0, it generates real and positive bound-states for the 
same value of a when q < (Fig. 1(b)). Further, Figs. 2(a) and (b) show the variation 
of the first three energy eigenstates as a function of a for (a) q = 1.0, V = 6m and (b) 
q = —1.0, V = 2m. Obviously, for the given V , as seen from Figs. 2(a) and (b) all possible 
eigenstates have negative (positive) eigenenergies if the parameter q is positive (negative). 
It is almost notable that there are some crossing points of the relativistic energy eigenvalues 
for some Vo values. 

The corresponding radial wave function ip nq (s) for the s-wave could be determined as 

Ms) = N nq s ib (l - qsfP^Xl - 2s), (47) 

where s = (1 + qe~ iax ) v . 

For the sake of comparing the relativistic and non- relativistic binding energies, we need 
to solve the ID Schrodinger equation for the complex form of the generalized WS potential 
given by Eq.(45). We employ the convenient transformation s(x) = (1 + ge -iax ) -1 , < r < 
00 — > < s < 1, to obtain 



1-2* „ ( HPs+jP-P 
(^)^ (S)+ (s-s 2 Y 



€ q (s) + 77— Ms) + 1 \ ^ ^ nq (s) = 0, (48) 
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for which 

f(s) = l-2s, a(s) = s-s 2 , or(s) = -(3 2 s + (3 2 - e 2 , 

e 2 = |^K (K<0), ^ = --1^^ (/? 2 >0). (49) 
n cr a erg 

Moreover, it could be obtained 



r(s) = -2(1 + c + e)s + (1 + 2c), c = ^/e 2 - /5 2 , (50) 

if 7r(s) = — (c + e)s + c is chosen for A;_ = — (c + e) 2 . We also find the eigenvalues 

A = - (c + e) (c + e + 1), A n = 2n (c + e + 1) + n(n - 1). (51) 

Finally, setting A=A n and solving for e, then the energy eigenvalues of the system under 
consideration could be found as 



E n (V ,q,ia) 



*2 2 r , -. n2 



2m 



n + 1 7 



(n + 1) 



, 7 = t^, 0<n<oo. (52) 



On the other hand, the radial wave function in the present case becomes 

V>n,(s) = V(l - s) e Pi 2c ' 2e) (l - 2s), (53) 

with s(x) = (1 + e~* a:E )~ 1 and 7V n9 is a new normalization constant determine by 

1 = ^2 , r (n + 26)!r(n + 2c+l) 2 f " (-l)y-* 

n<A J r(n + 2c + 2e + l) I £- Q p!(n - p)\(2e + p)\T(n + 2c - p + 1) 



x 



y (-i)vr(n + 2 C + r + 2e + i) | } s „ +2c+r _ _ +2edg 

2- r! (n-r)!r(2c + r + l) ji S U gSj dS ' ^ 



where 

i 



j s n+2c+r- P ^ _ qs yp+2e dg =2 F ^ n + 2c + r - p + i ^ p - 2e : n + 2c + r - p + 2 : l)B(n + 2c + r - p + 1, 





(55) 
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2. Non-Hermitian non-VT —Symmetric Generalized Woods- Saxon Potential 



Let two parameters; namely, Vo and q be complex parameters (i.e., Vo — > iVo, q — » ig), 
then we obtain the potential as 



V g (x) = Vo 



2cosh 2 (a:r) — sinh(2o;a;) — 1 


— i [cosh(a:;r) - 


sinh(ax)] 


1 + q 2 


2 cosh 2 (ax) - 


- sinh(2a::r) — 1 





(56) 



This potential is a non PT-symmetric but non-Hermitian possesses exact real spectra 
E nq = -^ + ^q 2 a 2 - 4V 2 - qa(2n + 1)) 



AV 2 + (jq 2 a 2 - AV 2 - qa(2n + 1)J 



16g 2 ' 
(57) 



if and only if 16g 2 m 2 > AV 2 + (Jq 2 a 2 -AV 2 - qa(2n + 1) j which gives - + + 

f (2n+l) < ^q 2 a 2 - 4V 2 . 

On the other hand, the corresponding radial wave functions i/) nq (s) for the s-wave could 
be determined as 

lM«) = N nq s b (l - S yP^\l - 2s), (58) 

i 

where s = (1 + ige _ax ) _1 . The integral I„ q (p,r) = J s n+2b+r ' p (l - s) p+2e ds is given by 

o 

I nq (p, r) = 2 Fi(n + 2b + r - p + 1, -p - 2e : n + 2b + r - p + 2; i)B(n + 2b + r - p + 1, 1). 

(59) 



3. Pseudo-Hermiticity and VT —Symmetric Generalized Woods-Saxon Potential 

When all the parameters V , a and q are complex parameters (i.e., V — > iV , a — > ia, 
we obtain 

^ = ~ g 2 + 2gshi (ax) + l [<? + + ' COS ^ = ~ x )- ( 60 ) 

This potential is a pseudo-Hermitian potential [27,48] having a 7r/2 phase difference with 
respect to the potential (I), rj = P-pseudo-Hermitian (i.e., PTV q (x)(PT)^ 1 = V q (x), with 
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P — 7] : x — > — x and T : i — > — i), it is also a non PT-symmetric but non-Hermitian 
having exact real spectrum given by 



Vr 



Enq = ~Tq + { ^ q2a " + + qai2n + 1} 



m 2 



16q 2 



4V 2 - ^q 2 a 2 + W 2 + qa(2n + 1)) 



2 ' 



(61) 

if and only if 16q 2 m 2 < 4V 2 - (^q 2 a 2 + 4V£ + qa{2n + 1)^ which gives -^0T)~2(^i) ~ 
f (2n+ 1) > ^q 2 a 2 + AV 2 

Figs. 3(a) and (b) show the variation of the ground-state (i.e., n = 0) as a function 
of the coupling constant Vb for different positive g with a = 10 and negative g with ct = 
I. Obviously, in Fig. 3(a), the P-pseudo-Hermitian non-PT-symmetric generalized WS 
potential generates real and negative bound-states for q > 0, it also generates real and 
positive bound-states for the same value of a when q < (Fig. 3(b)). Further, Figs. 4(a) 
and (b) show the variation of the first three energy eigenstates as a function of a for (a) 
q = 1.0, Vo = 2m and (b) q = — 1.0, Vo = 4m. Obviously, for the given Vo, as seen from 
Figs. 4 (a) and (b) all possible eigenstates have negative (positive) eigenenergies if the shape 
parameter q is positive (negative). There are some limitation on the number of bound states. 
It is seen that the range parameter a has limits for any given state. 

On the other hand, the corresponding radial wave functions ipnq(s) for the s-wave could 
be determined as 

Ms) = N nq s ib (l - S fPf b ' 2 -)(l - 2s), (62) 

i 

with s = (I + iqe-' iax )- 1 . The integral I ng (p,r) = J s n+2lb+r - p (l - s) p+2lt ds is given by 

o 

Inq(p, r ) =2 Fi(n + lib + r — p + 1, — p — 2ie : n + 2ib + r — p + 2; i)B{n + 2ib + r — p + 1, 1) 

(63) 



15 



IV. RESULTS AND CONCLUSIONS 



We have seen that the s-wave KG equation for the generalized WS potential can be solved 
exactly. The relativistic bound-state energy spectrum and the corresponding wave functions 
for the generalized WS potential have been obtained by the NU method. Some interesting 
results including the PT-symmetric and pseudo-Hermitian versions of the generalized WS 
potential have also been discussed for the real bound-states. In addition, we have discussed 
the relation between the non-relativistic and relativistic solutions and the possibility of 
existence of bound states for complex parameters. We have also shown the possibility to 
obtain relativistic bound-states of complex quantum mechanical formulations. Finally, the 
relativistic model provides real solution for the complex exponential potential where this 
solution is not available in the nonrelativistic model. 
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FIGURES 

FIG. 1. The variation of the ground-state (n = 0) energy with a = 1, in a non-Hermitian 
PT-symmetric potential given by Eq.(45), as a function of Vb for three different (a) positive and 
(b) negative q. 

FIG. 2. The variation of the first three energy eigenstates, in a non-Hermitian "PT-symmetric 
potential given by Eq.(45), as a function of a for (a) q = 1.0, Vq = 6m and (b) q = —1.0, Vq = 2m. 

FIG. 3. The variation of the ground-state (n = 0) energy, in a P-pseudo Hermitian 
non-PT-symmetric potential given by Eq.(60), as a function of Vb for three different (a) posi- 
tive q with a = 10 and (b) negative q with a = 1. 

FIG. 4. The variation of the first three energy eigenstates, in a P-pseudo -Hermitian 
non-PT-symmetric potential given by Eq.(60), as a function of a for (a) q = 1.0, Vq = 2m and (b) 
q = -1.0, Vb = 4m. 
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